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ABSTRACT 



Sensitive high angular and linear resolution radio images of the 240-pc radio 
jet in NGC 4151, imaged at linear resolutions of 0.3 to 2.6 pc using the VLBA and 
phased VLA at A21 cm, are presented and reveal for the first time a faint, highly 
coUimated jet (diameter ^1.4 pc) underlying discrete components, seen in lower 
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resolution MERLIN and VLA images, that appear to be shock-like features asso- 
ciated with changes in direction as the jet interacts with small gas clouds within 
the central ~100 pc of the galaxy. In addition, A21-cm spectral line imaging of the 
neutral hydrogen in the nuclear region reveals the spatial location, distribution 
and kinematics of the neutral gas detected previously in a lower resolution MER- 
LIN study. Neutral hydrogen absorption is detected against component C4W 
(E+F) as predicted by Mundell et al, but the absorption, extending over 3 pc, is 
spatially and kinematically complex on sub-parsec scales, suggesting the presence 
of small, dense gas clouds with a wide range of velocities and column densities. 
The main absorption component matches that detected in the MERLIN study, 
close to the systemic velocity (998 km s~^) of the galaxy, and is consistent with 
absorption through a clumpy neutral gas layer in the putative obscuring torus, 
with higher velocity blue- and red-shifted systems with narrow linewidths also 
detected across E+F. In this region, average column densities are high, lying in 
the range 2.7 x 10^^ Ts < Nn < 1.7 x 10^° Tg cm^^ K^^ (Tg is the spin tem- 
perature), with average radial velocities in the range 920 < Vr < 1050 km s~^. 
The spatial location and distribution of the absorbing gas across component E-l-F 
rules out component E as the location of the AGN (as suggested by Ulvestad et 
al.) and, in combination with the well-collimatcd continuum structures seen in 
component D, suggests that component D (possibly subcomponent D3) is the 
most likely location for the AGN. We suggest that components C and E are 
shocks produced in the jet as the plasma encounters, and is deviated by, dense 
clouds with diameters smaller than ~1.4 pc. 

Gomparison of the radio jet structure and the distribution and kinematics of 
ionized gas in the narrow hue region (NLR) suggests that shock-excitation by 
passage of the radio jet is not the dominant excitation mechanism for the NLR. 
We therefore favour nuclear photoionization to explain the structure of the NLR, 
although it is interesting to note that a small number of clouds with low velocity 
and high velocity dispersion are seen to bound the jet, particularly at positions of 
jet direction changes, suggesting that some NLR clouds arc responsible for bend- 
ing the jet. Alternatively, compression by a cocoon around the radio jet due to 
pressure stratification in the jet bow shock could explain the bright, compressed 
optical line-emitting clouds surrounding the cloud-free channel of the radio jet, 
as modelled by Steffen et al. 

Subject headings: galaxies: individual (NGG 4151) — galaxies: Seyfert — galax- 
ies: jets — radio lines: galaxies 
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1. Introduction 

Nuclear activity in galaxies is present over a wide range of luminosities, from the most 
distant and powerful quasars, to the weaker Active Galactic Nuclei (AGN) seen in nearby 
galaxies such as Seyferts, LINERs and even the Milky Way (e.g., Huchra & Burg 1992; Ho 
Filippenko & Sargent 1997). The standard model of nuclear activity involves the release 
of gravitational potential energy from galactic material accreted by a supermassive black 
hole at the galaxy center and, although radiation from the central engine is detected across 
the electromagnetic spectrum, it was the radio emission that first led to the discovery of 
powerful AGN (e.g., Baum & Minkowski 1960; Hazard, Mackey & Shimmins 1963; Bridle 
& Perley 1984). Indeed, the presence of powerful, highly-collimated relativistic radio jets in 
radio-loud quasars and radio galaxies, extending well beyond their host galaxy (e.g., Fanaroff 
& Riley 1974; Bridle & Perley 1984), provided key evidence for the exhaust material from a 
black-hole driven central engine (Scheuer 1974, Blandford and Rees 1974; Bridle et al. 1994; 
Urry & Padovani 1995). 

In contrast, radio-quiet quasars and Seyfert galaxies arc ten times more common but 100 
to 1000 times weaker at radio wavelengths than their radio-loud cousins (e.g., Goldschmidt et 
al. 1999); consequently nuclear starbursts have been advocated as the primary power source 
instead of black- hole accretion (e.g., Fernandes & Terlevich 1995). However, instrumental 
improvements and high resolution radio imaging with the VLA, MERLIN and the VLBA 
over the last twenty years have shown increasing evidence for coUimated radio emission 
in the form of small-scale radio jets that are weak analogues to jets in radio-loud AGN 
(e.g. Ulvestad & Wilson 1989; Wilson 1991; Nagar et al. 1999; Kukula et al. 1999) and 
indicate the presence of a central black hole and accretion disk, at least in some Seyferts. 
Similarly, measurements of high brightness temperatures in milliarcsecond resolution images 
of Seyfert nuclei with flat radio spectra have also suggested support for accretion-powered 
central engines (Mundell et al. 2000). 

Radio jets, although interesting in their own right, also provide a valuable probe of the 
interstellar medium close to the central engine, in particular the putative obscuring torus, 
which is thought to surround the central engine, providing fuel for the AGN and determining 
the observed optical spectral differences between type 1 (unobscured, viewed pole-on down 
the torus axis ) and obscured type 2 (viewed edge-on in the torus equatorial plane) (e.g., 
Antonucci & Miller 1985) . Some of the molecular gas in the torus is expected to be dissociated 
and ionized by the central UV/X-ray continuum source resulting in detectable quantities of 
neutral and ionized gas (e.g., Krolik & Lepp 1989; Pier & Voit 1995). Absorption of the 
radio jet continuum emission by intervening torus gas can be measured to determine ionized 
gas densities via free-free absorption (e.g., Ulvestad, Wrobel & Carilli 1999; Ulvestad et al. 



-4- 



1999), neutral gas columns via A21-cm neutral hydrogen (Hi) absorption (e.g., Mundell et 
al. 1995; Gallimore et al. 1999; Peck & Taylor 2002) and molecular gas content via OH 
absorption (e.g., Hagiwara et al. 2000). Due to the small size scales of the radio jet and the 
obscuring torus, high angular resolution imaging using the VLB A is required to spatially 
resolve gas in the absorbing region. 

In this context, NGC 4151, with its large quantities of neutral gas and elongated radio 
continuum structure, is an ideal candidate for such a study (see Ulrich 2000 for a compre- 
hensive review of the properties of the galaxy and its AGN). NGC 4151 is a Scyfcrt type 1.5 
(e.g., Ostcrbrock & Koski 1976) in an almost face-on {i=21°), grand-design, weakly-barred 
spiral host galaxy, Hubble type (R')SAB(rs)ab (de Vaucouleurs et al. 1991), that contains 
significant amounts of Hi throughout its two optically-faint spiral arms, fat weak bar and 
nuclear region (Davies 1973; Bosma, Ekers & Lequeux 1977; Pedlar et al. 1992; Mundell et 
al. 1995; Mundell & Shone 1999; Mundell et al. 1999). 

Radio continuum observations of the strong radio continuum nucleus show a linear radio 
structure in the form of a string of knots, elongated over ~3'/5 (230 pc) in average position 
angle (P.A.) ~77° (Wilson & Ulvestad 1982; Johnston et al. 1982; Carral, Turner & Ho 1990; 
Pedlar et al. 1993; Mundell et al. 1995), which is embedded in diffuse emission extending over 
~10'.'5 (~680 pc) (Johnston et al. 1982; Pedlar ct al. 1993). Carral et al. (1990) identified 
five main knots in the 3.5-arcsec jet at 15 GHz, annotating them CI to C5. A21-cm MERLIN 
observations of this central 4" region (Mundell et al. 1995) revealed localized and marginally- 
resolved Hi absorption, with a peak column density of A^h ~6x10^^ cm~^, against the 
component (C4) in the radio jet which is thought to contain the AGN; no absorption was 
detected against the other jet components to a hmiting column density of Ah ~2x10^^ Ts 
cm~^. An east- west column density gradient was observed and, in combination with UV 
column densities and early VLBI images, which showed C4 to consist of two components 
(C4E and C4W) separated by ~7 pc, led Mundell et al. (1995) to suggest that the weaker, 
western component, C4W, contains the optical/UV nucleus and thus the Hi absorption is 
taking place against the first component of the counterjct (C4E), due to gas in the obscuring 
torus. Structural and spectral index information obtained from subsequent radio continuum 
VLBA observations of the jet at 1.6 and 5 GHz led Ulvestad et al. (1998) to suggest a similar 
model but with the AGN located instead at the emission peak of C4E (component E after 
Ulvestad et al. 1998); requiring an absence of Hi absorption against the AGN, this model 
predicted Hi absorption only against the start of the counterjet, i.e. the tail of emission 
extending eastwards of the peak in C4E (component F after Ulvestad et al. 1998). 

In this paper we re-open the discussion on the possible location of the AGN based on new 
spectral fine observations, but whatever the precise location of AGN, NGC 4151 appears to 
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be unusual in having neutral gas located close to the nucleus, unlike other Seyferts that show 
H I absorption associated with gas located in dust disks on scales of 100—200 pc aligned with 
the host galaxy disk (e.g., Cole et al. 1998; Gallimore et al. 1999); in this respect NGC 4151 
seems more similar to Compact Symmetric Objects (CSO) and Steep Spectrum Core (SSC) 
objects (e.g., Conway 1996; Peck 1999; Peck & Taylor 2001). 

We present sensitive new, high resolution A21-cm continuum imaging, using the VLBA 
and phased VLA, of the entire 3'.'5 radio jet in NGC 4151 at an angular resolution of 40 mas 
(2.6 pc), along with higher resolution images of individual components with angular resolu- 
tions down to 5 mas (0.3 pc). In addition we present analysis of the first spatially-resolved 
images of the associated Hi absorbing gas with an angular resolution of 10.5 mas (0.7 pc) 
and confirm that the neutral gas, although clumpy, is located across the whole of component 
C4W (E-l-F after Ulvestad et al. 1998). Section 2 details the observational parameters and 
data reduction performed; in Section 3 we present the results of both the sensitive, high 
resolution continuum imaging of the jet and the H I distribution and kinematics. In Section 
4 we present the arguments for the location of the AGN being associated with component 
D and discuss the relationship between the radio jet and the narrow line region. We dis- 
cuss the Hi absorption in comparison with Lya absorption measurements and present our 
conclusions in Section 5. 

As discussed in Mundcll et al. (1999) the heliocentric radial velocity of NGC 4151 
is ~1000 km s^^, but distance estimates vary depending on the assumed value of Hq and 
whether the Virgocentric correction and the relative velocity of the Local Group with re- 
spect to the Virgo cluster are taken into account. Distance estimates lie in the range 
10<D<30.5 Mpc for 50<Ho<100 km s"^ kpc"^ and 1000<V<1523 km s'^, but since the 
uncertainty in the value of Hq is as large as the Virgocentric correction, for the purposes of 
this paper we assume a distance to NGC 4151 of 13.3 Mpc, using Ho=75 km s~^ kpc~^ and 
the heliocentric velocity of 998km s~^, giving a hnear scale of 0.065 pc mas~^ in the galaxy. 

2. Observations and Data Reduction 

The observations were obtained with the 10 element Very Long Baseline Array (VLBA 
- Napier et al., 1994) plus the 27 antennas of the Very Large Array (VLA - Thompson et al. 
1980) used in phased-array mode at 1.4 GHz, during a 14-hour observing run on 19 March 
1996. Phased-array mode involves the analog sums of the sampled and delayed intermediate 
frequency (IF) signals from the VLA antennas along each arm; the combination of the 
amplified and filtered signals from the three VLA arms is then recorded by the VLBI data 
acqusition system in VLBA format. When used in phased-array mode to augment the VLBA, 
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the VLA offers tfie equivalent sensitivity of a single 115-m antenna and an improvement 
in image sensitivity by a factor of ~2.4 compared with the VLBA alone. Although dual 
circular polarizations (right and left) were recorded, only the parallel hands (RR and LL) 
were correlated. The correlated data, using two bit sampling, consisted of two intermediate 
frequencies (IF) each with an 8 MHz bandwidth, 2 polarizations and 512 channels, resulting 
in a maximum spectral resolution of 15.6 kHz (3.3 km s~^) per channel. The first IF was 
centered near the redshifted Hi line of NGC 4151 at 1415.19 MHz corresponding to the 
recession velocity 1100 km s~^ (radio convention with respect to the Local Standard of Rest) 
while the second IF was shifted by 8 MHz relative to the center of IFl (center at 1423.19 
MHz) to provide a sensitive measure of continuum emission. The simultaneous VLA data, 
obtained with the array in C-configuration, were observed with a 25 MHz bandwidth centered 
at 1419.19 MHz encompassing the total 16-MHz VLBA bandwidth. 

Data editing and calibration followed standard methods^ and used the NRAO^ Astro- 
nomical Image Processing Software (van Moorsel, Kemball & Greisen 1996). VLBA ampli- 
tude scales were determined from standard VLBA antenna gain tables, maintained by NRAO 
staff, and measurements of Tsys made throughout the run; the integrated total flux density of 
NGC 4151, as measured by the VLA, was determined in the standard way from observations 
of 3C286. Observations of J1310-t-32 (VLBA antennas) were used to provide a check for 
relative amplitude differences between antennas and IFs, while J1642-f39 provided manual 
pulse and bandpass calibrations. Full phase referencing was not performed thus allowing 
the maximum possible observing time to be spent on NGC 4151, which is bright enough for 
self-cahbration. However short, regular observations of J1209-I-41 were interleaved through 
NGC 4151 scans with a cycle time of (11-1-2) minutes to provide a check on the delay-rate 
corrections; after imaging and phase/amplitude self-calibration, the amplitude-only correc- 
tions from J1209+41 were applied to NGC 4151. A single-channel pseudo-continuum dataset 
for NGC 4151, which was produced by averaging together the central 75% of both IFs, was 
then imaged and used as the input starting model for subsequent cycles of self-calibration 
and data editing (e.g. Walker 1995). The (m, f ) range covered by the data is large, 70 
kA — 40.8 MA at A21 cm, so the data are sensitive to a wide range of spatial frequencies. 
Using a combination of different data weightings and (it, v) tapers, it is therefore possible 
to produce a variety of images at different angular resolutions in order to show fully the 
complex structures present in the jet (Figures 1, 2, 3). Continuum images were formed with 



^Scc E. W. Greisen & P. P. Murphy 1998, The AIPS Cookbook, on-line at 

http: / /www. cv.nrao.edu/aips/cook. html. 

^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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angular resolutions from 40 mas (2.6 pc) to 5 mas (0.33 pc) and are discussed more fully 
in Section 3. Due to the lack of full phase referencing, absolute positions from J1209+41 
were not transferred to NGC 4151 and the R.A. and Dec axes are therefore provided by the 
pointing centre assumed correlation position, derived from earlier MERLIN observations; 
however, measured positions for E and D that are very close to those derived from the fully 
phase-referenced continuum images of Ulvestad et al. (1998), with the position of the peak 
flux density in component E measured, from our 5-mas uniformly weighted continuum image 
(Figure 3), to be a = 12'^10"*32^5821, S = 39°24'21'.'060 compared with a = 12'^10"*3215822, 
5 = 39°24'21'.'059 from Ulvestad et al. (1998). 

The additional data editing and corrections derived from the self-calibration of the 
continuum emission of NGC 4151 were finally applied to the NGC 4151 spectral line data 
and all subsequent processing concentrated on IFl. The continuum contribution in IFl 
(formed from the line-free channels 32-182 and 350-450) was subtracted in the (w, v) plane 
using AlPS task UVLIN and the corresponding 'pure continuum' dataset was formed from the 
same line-free channels in the {u, v) plane using AVSPC. Concentrating on the region of the 
source where absorption is expected (i.e., components D, E, F), the pure continuum and 
line datasets extracted from IF 1 were then Fourier transformed with uniform and natural 
weighting (with 1-mas pixels), deconvolved and added together to form the corresponding 
absorption cube (256 x 256 x 512) with the correct continuum level adjacent to the fine. It 
should be noted that although a combination of data from IFl and IF2 were used to produce 
maximum-sensitivity continuum images of the whole jet, this continuum image is not suitable 
for addition to the spectral line dataset (from IFl) to form optical depth maps because small 
errors in the bandpass calibration (<5%) between the two IFs produce significant errors in 
optical depth spectra (Faison et al. 1998); therefore only continuum emission from line-free 
channels in IFl were used to produce the optical depth maps. A full spectral-resolution cube 
(spectral resolution 3.3 km s~^) with 512 channels was produced, whilst spectral averaging 
was used to increase the channel sensitivity in spectral line cubes imaged with 128 and 256 
channels. The intermediate spectral-resolution natural cube with 256 channels (6.6 km s~^ 
per channel) was used as a consistency check to confirm the presence of multiple narrow 
components detected in the 128- and 512-channel cubes (Figure 4). The more sensitive 128- 
channcl cube (spectral resolution 13.2 km s""*") was imaged and used to construct optical 
depth maps (see Figure 5). which were created using only absorption stronger that 1.5 mJy 
beam~^ (3o"). In Figure 5, structures much smaller than the beam are an artefact of the 
data processing and are compatible with the signal-to-noise of the spectra; trends on scales 
larger than or comparable to the beamsize are real. Spectra extracted from a higher angular 
resolution, uniformly-weighted cube indicates clumping of the gas on scales ~6 mas, but the 
signal-to-noise of the spectral was insufficient for reliable moment analysis. 
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3. Results 

3.1. A21-cm Radio Continuum Emission 

As discussed in Section 2, the A21-cm (1.4-GHz) VLBA+VLA data are sensitive to a 
wide range of spatial scales, allowing a variety of images at different angular resolutions to 
be produced, fully probing the structures in the jet. Figure 1 shows continuum emission 
from the entire 3'.'2-long jet (~210 pc) imaged at 40 mas (2.6 pc) resolution. For the first 
time, weak emission from a highly-collimated radio jet is seen underlying the five discrete 
components, or knots (CI— C5), seen in previous lower resolution MERLIN and VLA images 
(e.g., Carral et al. 1990; Pedlar et al. 1993). These five knots are now spatially resolved and 
can be clearly seen, in the the 25-mas resolution image (Figure 3 - upper panel), to coincide 
with changes in direction of the jet. 

Due to the large number of features visible along the jet in the VLBl images, the 
"CI— C5" labeling scheme used for MERLIN and VLA images, as noted by Ulvestad et al. 
(1998), is no longer adequate; we therefore adopt the nomenclature introduced by Ulvestad 
et al. (1998) in which the components are labeled from west to east with letters from A to 
H. We identify two new additional components, I (east of H) and X (between C and D). 
Table 1 lists the correspondence between the component designations used in Mundell et al. 
(1995) and Ulvestad et al. (1998). 

Connecting components I and A with a single straight line gives an average jet P.A. of 
77°, as found in lower resolution studies (e.g. Pedlar et al. 1993). However, the jet structure 
is more complex than suggested by a single P.A and, as can be seen in Figures 1 and 3(a), 
several sections of the jet point in a predominantly east-west direction, along P.A. ~89°, i.e. 
sections connecting components A to the lower tip of B, F to the lower tip of G and the 
eastern end of H to 1; finally and most obviously, the peaks of the three subcomponents in 
D (Figures 2, 3(b)) lie along a line in P.A. ~84° and connect with X. However, these almost 
horizontal sections are not colinear, as the jet exhibits abrupt and almost step-like, direction 
changes. These changes are most obvious in Figure 1 at components B, G and H, and in 
Figures 2 and 3(b) at component E. The subarcsecond-scale jet structures therefore do not 
align with the average P.A. 77° and it is clear from Figure 1 that the final, 'emergent' P.A. 
for both eastern and western ends of the jet after it has finished its various deviations is 
indeed in P.A. ~90°. 

Although the exact location of the AGN between F and C is not known (see Section 4 for 
discussion), the jet appears initially very well-coUimated in the vicinity of the components 
associated with D (Figure 2). Despite deviating quite abruptly at E, the eastern jet (through 
components F, G, H, I) remains more highly collimated along its length than the western 
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jet (through components C, B, A), which contains more diffuse and extended components. 
No emission is detected in the region between component X and the eastern tail of C; to a 
3-0" detection threshold of 0.15 mJy beam~^, which corresponds to brightness temperature 
upper limit of Tb ^ 2.8x10^ K at A21 cm in the 40-mas resolution image. 

3.2. Parsec-Scale Neutral Hydrogen Absorption 

Neutral hydrogen absorption is only detected against continuum components E+F; no 
Hi absorption is detected against component D to a limiting 3-a absorption depth of 1.5 
mJy beam~^, i.e. a maximum column density A^h ~1.1x10^^ cm~^ (per 13.2 km s~^ 
channel), assuming a constant column density across D (sec Section 4 for discussion). Fig- 
ure 4 shows some representative absorption spectra taken at five locations (labelled a to 
e) across the 3-pc-long continuum structure of the 'banana-shaped' component, E-l-F; the 
spectra were extracted from the 512-channel naturally- weighted cube with angular resolution 
10.5 mas (0.7 pc) and spectral resolution 3.3 km s~^. For the purposes of moment analy- 
sis, we derived the distribution of column density and corresponding gas kinematics from a 
spectrally-averaged, 128-channel (spectral resolution 13.2 km s~^), naturally-weighted cube 
with a restoring beamsize 11.5x9.0 mas (0.75x0.59 pc) in P.A. 0°. Only absorption stronger 
than the 3-cr absorption depth of 1.5 mJy beam~^ per channel was included in the moment 
analysis. 

3.2.1. Column Densities, Gas Distribution and Kinematics 

The distribution of H I column density (Ah) per unit spin temperature across compo- 
nents E+F is shown in Figure 5(a). As can be seen in Figures 4 and 5, Hi absorption is 
detected against the whole extent of E+F (C4W) and is spatially and kinematically complex; 
across E+F, average column densities lie in the range 2.7 x 10^^ Ts < Ah < 1.7 x 10^*^ Tg cm~^ 
where Ts is the spin temperature for which the value is unknown but is typically 
~10^— 10^ K in the Galaxy (Heiles & Kulkarni, 1988) and average radial velocities lie in 
the range 920 < Vr < 1050 km s~^. These values are consistent with the lower resolu- 
tion (0'.'16x0'.'15x26.3 km s"^) MERLIN absorption study in which Mundell et al. (1995) 
reported an E-W column density gradient, across the marginally resolved component C4, 
3 X 10^^ Ts < Ah < 6 X 10^^ Ts cm'^ K'^ and an absorption line with Full Width at 
Half Maximum (FWHM) 91 km s~^ centered at 998 km s~^. The higher column density 
values measured in the VLBA datacubc confirms that the absorption is spatially resolved 
and occurs on scales significantly smaller than the MERLIN beam. 
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The mean column density gradient decreases towards the NE, along the tail of F, but 
higher column- density clumps arc present throughout this region, as measured from a higher 
resolution uniformly- weighted cube (not shown); in particular, the peak column density A^h 
~1.7 xlO^° Ts cm~^ does not occur at the location of the continuum emission peak (E) but 
instead occurs ~15 mas (1 pc) NE of the southern tip of the 'banana' (or ~0.5 pc from the 
continuum emission peak). Column densities in the region of the continuum peak are still 
significant, but lower at 9 x lO^'^ Tg < A^h < l-2x lO^o Tg cm-^ (Figure 5(c)). 

Figure 5(b) shows the mean velocity field derived from single Gaussian moment analysis. 
Although useful to show the general trend in velocity across the 'banana', this velocity field 
represents only the main velocity component close to the galaxy systemic value. As can 
be seen in Figure 4, multiple, narrow velocity components are present at different positions 
across the continuum structure. At this high spectral and spatial resolution one can see 
a family of complex absorption lines (e.g. location a), multiple velocity components (e.g. 
location /?) and narrow single lines (e.g. locations 5, e) at different locations across the 
continuum structure. The velocity components are indicated in Figure 4 and their properties 
are listed in Table 2. 

4. Discussion 

4.1. The Location of the AGN 

The location of the AGN in NGC 4151 is broadly accepted to be in component C4 
(e.g. Pedlar et al. 1993; Mundell et al. 1995). However, as discussed previously, it is now 
known that C4 consists of a number of subcomponents and the precise association of the 
AGN with one of these subcomponents is still a matter of debate. In this study, we use the 
characteristics of the radio continuum structure and the H I absorption to argue for the AGN 
located in component D. 

4.1.1. Radio Continuum Structures 

Components C and (E+F) share some structural similarities having a 'banana'-like ridge 
on their northern half, but with the tail of (E-l-F) pointing eastwards (i.e. counterjet) and 
that of C pointing westwards (forward jet); if these opposed arc-like structures are due to 
interactions of the radio jet with dense clouds in the circumnuclear environment, making 
each component an impact site, their similar shapes and opposing directions suggest that 
the AGN hes somewhere between the two. The greater flux density of E-l-F suggests that 
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the AGN is located closer E+F than to C, or that C is an older component. In addition, 
as described in Section 3, the continuum structure of component D shows a high degree of 
collimation; indeed the structure resembles radio core-jet structures seen in higher power 
AGN such as radio loud quasars. If the AGN corresponds to one of the components in D 
(possibly D3), as originally suggested by Mundcll et al. (1995), this would explain the very 
linear, core-jet-like, structure of D and the greater brightness of E+F relative to C. 

Assuming that D delineates the well-collimatcd jet before it has suffered any interference, 
a jet diameter of approximately 1.4 pc is inferred. A small dense gas cloud that wanders 
into the path of the radio plasma flow would disrupt or deviate the jet flow. The amount of 
disruption depends on the impact parameter of the jet with the cloud, assuming the cloud 
diameter is less than that of the jet (<1.4 pc). Component E-l-F can then be understood 
as a grazing encounter with a small dense cloud, so the (counter-)jet is deflected at E but 
continues to flow over the cloud to F and outwards through G, and retaining some degree of 
collimation through to I. 

To the west of D, the similar but more diffuse, but clearly shell-like structure of compo- 
nent G (Figure 2) can be explained by the jet flow to the south-west of the AGN impacting 
a cloud more directly and so being more disrupted. With a cloud diameter smaller than the 
jet diameter, the jet flow is not terminated at C but continues to flow around the cloud and 
onwards to components B and A, but with signiflcantly more disruption and poorer collima- 
tion than its counterpart to the north-cast. Indeed, the flattening of the northern half of the 
shell structure of C might be explained by the containment effect of a dense bounding cloud 
to the north of C. As shown by HST studies of ionized gas clouds in the NLR of NGC 4151 
(e.g., Winge et al. 1997; Hutchings et al. 1998; Kaiser et al. 2000), a large number of gas 
clouds exist in the inner regions of the galaxy, so the chance of the radio jet intercepting 
and being affected by gas clouds is high. In addition, the relatively flat radio spectrum 
of component E might be explained by re- acceleration of synchrotron-emitting electrons in 
the shock produced by the interaction with the cloud, rather than the AGN-interpretation 
suggested by Ulvestad et al. (1998). 



4.1.2. Distribution of Hi Absorption 

Mundell et al. (1995) originally suggested that the discrepancy between the high column 
densities and systemic gas velocity inferred from MERLIN radio H I absorption measurements 
(peak A^H ~6.0 x 10^^ cm^^ at ~998 km s^^ for Ts=100 K ) and much lower columns inferred 
from UV observations of Lya absorption (Kriss et al. 1992) with signiflcant blueshifts 
(6 X lO^'' < Ah < 6 X 10^° cm~^) could be interpreted as support for the presence of 
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the putative obscuring torus invoked in Unification Schemes (e.g. Antonucci Sz Miller 1985). 
In this interpretation, the lower column densities of blueshifted gas seen in the UV are 
interpreted as outflowing ionised gas seen along the line of sight to the AGN; the higher 
column densities and lower gas velocity seen in the radio are interpreted as absorption 
through gas in the obscuring torus along the line of sight towards the first component of the 
counter] et. Component C4W (or D) was then assigned as the true optical/UV nucleus with 
C4E (E+F) representing the first radio component in the counterjet against which the H I 
absorption was occurring (Mundell et al. 1995). 

Ulvcstad ct al. (1998) used the VLBA to image, at high angular rcsohition, the bright 
components (C, D, E, F) at A18 cm and A6 cm and set upper limits for the emission at 
A2 cm. They revised the model of Mundell et al. by suggesting an alternative location 
for the AGN in E (as discussed above) and predicted that Hi absorption would only be 
detected against subcomponents of F (F3— F5 in Ulvestad et al. 1998), with subcomponents 
E— F2 lying behind fully ionized gas and therefore resulting in no Hi absorption against 
these components. This prediction is easily testable with the present data and found to be 
inconstent; Hi absorption is detected across the whole extent of E+F (see Section 3 and 
Figures 3, 4), ruling E out as the location for the AGN. 

Considering both the Hi absorption results and the radio continuum structures (as 
discussed above), component D remains the most likely location for the AGN. Assuming 
an ionized/neutral/molccular 'onion-skin' model for the torus, component D would then lie 
behind the ionized inner skin of the torus with E+F lying behind the neutral (dissociated) 
gas layer and molecular gas would be further out. The geometry of such a model is illustrated 
in Figure 6. In this model, the weak diffuse component X, which is resolved out at higher 
resolutions and radio frequencies, has been placed at the centre of the torus since it also 
lies along the coUimation axis defined by D and so could be an alternative location for the 
AGN; the jet can be shifted relative to the H2 image (relative astronometry is poor) to 
place whichever radio component contains the AGN at the torus centre, with the condition 
that component E+F must lie behind the neutral gas layer. If subcomponent D3 contains 
the AGN and ionized gas fills the inner region of the torus between D3 and the neutral 
gas layer starting at E, the radius (Ri) of the ionized gas sphere interior to the neutral 
layer would be approximately Rj ~3.3(sin^^^^) pc (where 9 is the angle between the jet axis 
and the line of sight). For an ionizing fiux of 1.5x10^^ photons per second (Kriss et al. 
1995; Ulvestad et al. 1998), the electron density in Stromgren radius Ri ~3.3(sin~^^) pc 
is ne~3.8xl0^ (sin^/^^) /~°'^ cm~^ (where / is the fiUing factor); this electron density is 
intermediate between the densities expected for Seyfert narrow-line regions (~10^ cm~^) 
and broad-line regions (~10^ cm^'^) and similar to values inferred from free-free absorption 
in Mrk 231 and Mrk 348 (Ulvestad et al. 1999) and thermal bremsstrahlung emission upper 
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limits in NGC 4388 (Mundell et al. 2000). 

Sensitive radio imaging over a range of wavelengths to derive radio spectra of the sub- 
components in D would help to identify a fiat-spectrum, high brightness temperature com- 
ponent indicative of an AGN and characterise the quantity of ionized, free-free absorbing 
foreground gas. Alternatively, the nucleus itself may not be visible at centimetric radio wave- 
lengths if significant quantities of ionized gas lie along our line of sight and free-free absorb 
radio photons (e.g. the absence of A21-cm continuum emission between components X and 
C can be seen in Figure 1 and is described in Section 3); sensitive, high frequency radio 
observations, radio proper motion studies and improved HST and VLBA relative astrometry 
will be important in determining the true location of the AGN. 

4.2. The Relationship between the Ncirrow Line Region and the Radio Jet 

Ionized gas in the narrow line region of NGC 4151 is distributed in a clumpy, biconical 
structure, extending over ~2" north-east and south-west of the bright nucleus, and has been 
modelled as gas outflowing at a constant velocity of ~350 km s^^ along the surface of a 
bicone with total opening angle of 80°, viewed at 10° outside the cones (Hutchings et al. 
1998). The blueshifted gas southwest of the nucleus is approaching (see also Schultz 1990; 
Evans et al. 1993; Robinson et al. 1994; Boksenberg et al. 1995). 

A key question in the study of the energetics and kinematics of gas in NLR regions of 
Seyferts is the importance of the radio jet. Competing models to explain optical line emission 
from the NLR advocate two main excitation mechanisms; direct photoionisation by the 
central continuum source (e.g., Binette et al., 1996) or shock ionisation caused by the passage 
of the radio jet through the ISM (e.g., Dopita & Sutherland, 1995, 1996; Bicknell et al, 1998). 
Wedge-shaped morphologies observed in some NLR's (e.g., Evans et al., 1991; Wilson et al., 
1993; Wilson & Tsvetanov, 1994; Simpson et al., 1997) are consistent with photoionisation 
by a cone of UV radiation from the central engine. Alternatively, the correlation between 
radio power and [Om] luminosity and line-width (de Bruyn & Wilson 1978; Wilson & Willis 
1980; Whittle, 1985, 1992), and the approximate spatial coincidence between radio and 
optical structures in other NLRs suggest that the radio jets may be responsible for some of 
the optical line emission. 

Radio observations of powerful quasars and radio galaxies have shown evidence for bulk 
relativistic (non-thermal) flows in the radio jets (e.g. Pohl et al. 1995; Hough, Zensus & 
Porcas 1996 ), but it remains controversial whether ejected radio plasma in Seyfert galaxies 
has bulk relativistic motion (e.g. Ulvestad et al. 1999; Brunthaler et al. 2000). Bicknell et al. 
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(1997) suggest that Seyfert jets may be mildly relativistic close to the black hole and become 
increasingly dominated by thermal plasma on larger scales, as they entrain material from the 
ISM. They argue that the energetic importance of Seyfert radio jets has been underestimated 
and that jets may provide the energy input to drive the optical emission in the NLR, with 
jet energy fluxes (based on a study of NGC 1068) being greater than previously supposed. 
Nevertheless, a direct causal relationship between the radio jet and NLR emission has not 
yet been established and the relative contributions of direct photoionisation of emission hne 
regions by the central continuum source, and strong interaction of the emission-line gas with 
radio plasma remains controversial. 

Until now, the relatively small angular extent (~ few arcsecs) of Seyfert radio jets has 
limited the detail visible in VLA and MERLIN images and, radio images at even r-'OfH angular 
resolution, have proven insufficient to identify a clear association between NLR clouds and 
radio jet components (e.g. Mundell et al. 1995; Kaiser et al 2000). Now however, the 
improved sensitivity of new radio VLBI studies is revolutionising the study of Seyfert radio 
jets. Figure 7 shows our 40-mas resolution image of the radio jet in NGC 4151 (from Figure 
1) overlaid on the F501N narrow-band image taken with the Faint Object Camera (FOC) on 
board the Hubble Space Telescope (HST) (image courtesy of R. Catchpole; see also Winge et 
al. 1997). The FOC image with 0'.'0143 pixels has been resampled to match the 0'.'0175 pixel 
size of this radio image and, since relative astrometry between the HST and VLBA images is 
insufficient for accurate 'blind' alignment, the peak of the optical nuclear emission (assumed 
originate from the AGN) was aligned with the peak of radio emission in component D, 
assumed to contain the AGN. A rehable comparison of the projected location of the optical 
line-emitting clouds and radio jet features is now possible, although any physical association 
along the line of sight is still unknown. 

The misalignment of the radio jet and NLR supports photoionization of the NLR by 
the AGN as the dominant excitation mechanism. A similar misalignment is present between 
soft X-ray emission imaged with Chandra (Yang, Wilson & Ferruit 2001) and the radio jet. 
However, a small number of bright [OlIl]-emitting clouds appear to be closely associated with 
the jet, lying close to and appearing to bound the radio jet emission; these clouds are labelled 
6, 9, 11, 12, 15, 17, 19, 20 in Figure 7. Kaiser et al. (2000) classify all of these clouds as low 
velocity (1^1 < 400 km s~^) but high velocity dispersion {Av > 130 km s~-^) and suggest that 
they may be intrinsically high-velocity clouds with smaller projected components of velocity 
along the hne of sight. In addition. Kaiser et al. (2000) note that clouds 15, 19 and 20 also 
have high velocity, blueshifted clouds {\v\ > 1000 km s~^) associated with them that are too 
faint to be visible in their labelled WFPC2 images. It is therefore likely that the radio jet 
might have cleared a channel through the NLR and suffered disturbance as a consequence; 
remaining NLR clouds then act to gently deviate the plasma flow. Alternatively, as modelled 
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by Steffen et al. (1997), the enhanced optical emission from the clouds bounding the radio 
jet might arise due to compression by a cocoon around the radio jet; a structure of cylindrical 
nested shells are produced around the radio jet due to pressure stratification in the jet bow 
shock and the brightest compressed clouds then surround the cloud-free channel of the radio 
jet (see Figure 2 in Steffen et al. 1997). Numerical simulations, by Steffen et al., of this jet 
propagating through an inhomogeneous NLR-like medium show optical linewidths as high as 
500 km s~^, and large variations in the cloud radial velocities, qualitatively similar to those 
observed for NGC 4151. 



4.3. Complex Hi Absorption and the Obscuring Torus 

4-3.1. Comparison with UV/X-ray Absorption 

Absorption by gas along the hne of sight to AGN is observed at a variety of wavelengths, 
with X-ray and UV absorption occurring predominantly in ionized gas and radio absorption 
in neutral and molecular gas. The nature of the "warm" (or ionized) absorbing gas detected 
in X-ray spectra of the nuclei of nearby AGN is still unknown. As many as 50% of nearby 
AGN show absorption by ionized gas (e.g. George et al. 1998) with many also show associ- 
ated UV absorption by highly ionized species such as C IV (Crenshaw et al. 1999; Kraemer 
et al. 2001). In NGC 4151 broad, saturated Lya absorption lines are detected in the UV 
(Kriss et al. 1995; Espey et al. 1998) and, when considered along with the properties of 
strong absorption lines due to other species such as Nv CiV Si IV, indicate that the absorber 
does not fully cover the source of background emission (Kraemer et al. 2001). The ionized 
gas appears to be outflowing from the nucleus with radial velocities as high as —1680 km s~^ 
(Kraemer et al. 2001) and cannot account for the observed X-ray absorption (e.g. George 
et al. 1998; Yang et al. 2001). Kraemer et al. (2001) suggest a model in which the UV 
absorber covers the BLR, with the X-ray absorber closer to the central source and obscuring 
a smaller X-ray emitting region. 

The Hi column densities in NGC 4151 derived from Lya absorption (e.g. Kriss et al. 
1992, 1995; Espey et al. 1998), ~3 x 10i^<A^H(Lya) <6xl020 cm'^, are significantly lower 
than those derived from the present A21-cm observations (see Section 3), and although the 
UV continuum is known to be highly variable, the UV- measured column densities remain 
lower than the radio values. This conflrms that the neutral gas obscuring the radio continuum 
emission is distinct from the UV-absorbing gas and we re-propose our geometrical model 
(Figure 6) in which the A21-cm Hi absorbing gas is located in a neutral, photo-dissociated 
layer of the obscuring torus, detected in front of the flrst radio component in the radio 
counterjet (E-l-F), while the UV-absorbing gas is outflowing, ionized gas (possibly blown 
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from the torus walls) along the hne of sight to the nucleus. 

4.3.2. A Clumpy Torus? 

The total absorption spectrum across E+F is consistent with that detected in the MER- 
LIN study (Mundell et al. 1995), but the presence of additional narrow absorption lines with 
higher velocities (Figure 4, Table 2) across E+F argues for dense clumps of gas within the 
circumnuclear disk or torus, and may indicate turbulence within the gas. dumpiness of 
torus gas is predicted theoretically (e.g. Krolik & Bcgclman 1988), with cloud covering 
factors close to or greater than unity, and recent X-ray determinations of column density 
variations of 20%— 80% over timescales of months to years (e.g. Risaiti, Elvis & Nicastro 
2002) in Scyfcrt 2 nuclei support clumping of torus material on scales less than 1 pc. Clouds 
drifting across the background continuum source would then explain observed column den- 
sity variations. Sub-parsec scale radio monitoring of the H I absorption across component E 
in NGC 4151 is required to determine if the neutral gas columns vary, and might therefore 
be associated with drifting cloud cores, or remain constant, suggesting Hi it is the primary 
torus volume filler in which denser clouds are embedded. 



4.3.3. NGC 4151 as a Nearby Radio Galaxy? 

A similar scenario to that suggested for NGC 4151 has been proposed for Hi absorption 
detected in the nuclear regions of the compact symmetric object 1946-1-708 (Peck, Taylor & 
Conway 1999), in which complex Hi absorption is detected against the counterjet (see also 
Taylor 1996; Taylor et al 1999). Indeed the higher rate of detection of Hi absorption in 
AGN with two-sided, non-boosted radio jets such as compact symmetric objects (CSO) and 
compact steep spectrum objects (CSS) compared with core-dominated radio sources, with 
Doppler-dimmed counterjets, provides support for the Unification Scheme in which objects 
seem close to the plane of the sky have a more favourable geometry and orientation for the 
detection of HI absorption against the counterjet due to gas in a nuclear obscuring torus (e.g. 
van Gorkom et al. 1989; Pihlstrom 2001). NGC 4151 is an unusual Seyfert in having neutral 
gas located close to the nucleus, unlike other Seyferts that show H I absorption associated 
with gas located in dust disks on scales of 100—200 pc ahgned with the host galaxy disk 
(e.g.. Cole et al. 1998; Gallimorc et al. 1999); in this respect, the Hi absorption, and indeed 
the well-coUimated radio jet in NGC 4151, closely resemble parsec-scale structures in CSO 
and CSS. 
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5. Conclusions 

We have used the VLB A and phased VLA at A21 cm to study the parsec-scale radio 
continuum emission and distribution and kinematics of neutral hydrogen in absorption in 
the Seyfert galaxy NGC 4151. We find: 

• A faint, highly coUimated radio jet, diameter ^1.4 pc, underlying brighter radio knots, 
seen in lower resolution MERLIN and VLA images, that appear to be shock-hke fea- 
tures associated with changes in direction as the jet interacts with small gas clouds 
within the central ~100 pc of the galaxy. 

• H I absorption across the whole of component C4W (E+F), as predicted by Mundell 
et al. (1995), with complex spatial and kinematic structure. 

• The location and distribution of the HI absorption and the structure of the radio 
continuum emission argue for the AGN being located in component D, rather than 
component E as argued by Ulvestad et al. (1998). 

• We suggest the absorbing gas lies in a thin photo-dissociated layer of clumpy neutral 
gas in between the molecular and ionized gas in the circumnuclear torus, with an inner 
radius ~3.3 (sin~^^) pc. 

• The location of the absorbing gas close to the AGN is unusual for a Seyfert galaxy and 
we argue that the H I and continuum properties of NGC 4151 arc similar to those seen 
in compact symmetric objects and compact steep spectrum objects. 

• Spatial association of the radio jet knots and optical line-emitting clouds with high 
velocity dispersion in the narrow line region suggest some interaction between the radio 
jet and the clumpy interstellar medium, but photoionization remains the dominant 
excitation mechanism for the NLR. 
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Table 1: Radio Component Designations (see also Figure 1) 

Mundell et al. (1995) Ulvestad et al. 1998 

and This Paper 

B 

C 

xt 

D 
E 
F 
G 
H 

It 

^New component identified in this paper 



CI 

C2 

C3 

C3-C4, eastern bridge 

C4W 

C4E, peak 

C4E, eastern 'tail'. . . 
C5, western extension 

C5 

C5, eastern extension 
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Tablc 2: Central velocity, width and mean column density for absorption components de- 
tected at five locations across E+F, as shown in Figure 4. Systemic velocity is 998 km s~^. 



Box Label 


Component Velocity 


Width 


Mean A^h/Ts 


(see Fig. 1) 


((3ptical Heliocentric km s^ '" 


) (kms-i) 


(10^-" cm-2 K-i) 


LX 


Q98 7 


16 


7 


a 


962.2 


20 


1.1 


a 


1022.3 


26 


1.0 


a 


1132.5 


20 


0.6 


LX 


147fi 

J-^ / u.u 


oo 


1 9 




825.2 


20 


1.8 




925.4 


20 


1.8 




978.9 


20 


2.1 




1132.5 


20 


0.9 




1443.2 


20 


0.9 


7 


848.6 


13 


0.8 


7 


982.2 


20 


2.0 


7 


1443.2 


20 


1.2 


5 


1002.2 


40 


4.8 


5 


1323.0 


20 


1.4 


5 


1673.7 


20 


1.7 


e 


988.9 


26 


3.6 


e 


1045.7 


20 


1.1 


e 


1209.4 


20 


1.1 
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Fig. 2. — Naturally- weighted VLA+VLBA A21-cm continuum image showing components 
C, X, D, E, F at high angular resolution. Contour levels are (-0.25, 0.25, 0.5, 1, 2, 3, 4, 5) 
mJy beam~^ and the restoring beam, with size 12.7 x 9.0 mas (0.83 x 0.59 pc) in P.A. 16°, is 
shown in the lower left corner of the image. As in Figure 1, Ulvestad et al. (1998) component 
labeling convention is shown above the jet (plus new component X identified in this paper) 
and corresponding component labeling from lower-resolution studies (e.g., Harrison et al. 
1986; Mundell et al. 1995) is shown below the jet (C3, C4E, C4W). 
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Fig. 3. — Upper panel: full 3'.'2-lengtli jet imaged with restoring beamsize of 25 mas (1.6 pc). 
Lower panel: maximum angular resolution, uniformly weighted, image with restoring beam- 
size of 5 mas (0.33 pc - 1.1 light year), with components marked with Mundell et al. (1995) 
convention below and the Ulvestad et al. (1998) convention above. 
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Fig. 4. — Naturally- weighted 1.4-GHz continuum image with examples of spectra across this 
component; the spectra are taken from the naturally-weighted full spectral resolution (3.3 
km s~^) data cube and are displayed with 3-channel (FWHM) Gaussian smoothing applied. 
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(c) Nh (10^" cm -2) (d) Nh (10^" cm 2) 




Fig. 5. — Panels showing moment maps of Hi column density (A^h) and velocity along the 
'banana' component E+F: (a) A^h map with A^h contours superimposed; the peak A^h is 
1.7 X 10^° Ts cm~^(b) velocity map with velocity contours superimposed (c) A^h map with 
contours of 1.4-GHz radio continuum emission superimposed with contour levels (-0.8, 0.8, 
1.6, 3.2, 4.8, 6.4) mJy beam~^ (d) A^h contoured on velocity. The 11.5x9.0-mas (0.7x0.6 
pc) beam (P.A. 0°) is shown in the lower left corner of each image, but omitted from (c) for 
clarity. 
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Fig. 6. — Montage of the inner 250-pc of NGC 4151, showing torus of H2 emission in green 
(from Fernandez et ah 1998), ring of Hi inferred from absorption measurements in blue and 
1.4-GHz radio continuum emission from radio jet in red. Ionized gas (black) is assumed to 
fill the torus inside the H I ring. The segments missing to the north and south of the H2 torus 
are due to the limited filter width that excluded the high velocity wings of the H2 line; this 
provides evidence for rotation of the torus as the northern and southern segments contain 
the gas with the highest radial velocities (i.e. line wings) if the torus is rotating. 
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